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The Carbon Isotope Biogeochemistry of EC0 2 Production in a 
Methanogenic Marine Sediment 

ABSTRACT 

To investigate the relationship between EC0 2 fi^C values and rates of the 
dominant remine raliza tion processes at the organic-rich field site of Cape Lookout 
Bight, NC, the isotopic composition of porewater EC0 2 was measured on a seasonal 
basis. The EC0 2 S 13 C values varied seasonally in response to changes in rates of 
sulfate reduction and methanogenesis, the dominant remineralization processes at this 
site. 

A tube incubation experiment was also performed to determine the isotopic 
signature of the EC0 2 produced by sulfate reduction and methanogenesis. The 6 C 
of the EC0 2 produced in the sulfate reduction zone determined from the tube 
incubation was -14.3 ± 1.9, a value enriched in ^ 3 C relative to the labile organic 
fraction. The l3 C-enrichment may be caused by low rates of methanogenesis 
occurring in the sulfate reduction zone. The 6 13 C of the EC0 2 produced in the 
methanogenic zone was estimated to be +44 per mil, whereas the co-produced 
methane was -65 per mil. The fractionation factor for C0 2 reduction was calculated 
to be 1.055, a value in agreement with previous estimates at this site. The measured 
concentration and 6 13 C of the EC0 2 at Cape Lookout was closely reproduced by a 
diagenetic model using the measured rates of sulfate reduction and E C0 2 production, 
and the isotopic signature of the EC0 2 production in the two biogeochemical zones. 
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INTRODUCTION 


The isotopic composition of total dissolved inorganic carbon (EC0 2 ) in 
sediments is a unique record of diagenetic processes. Presley and Kaplan (1968) used 
the isotopic signatures of EC0 2 to confirm that the downcore increases in EC0 2 from 
nearshore sediments were a result of metabolic activity of organisms in the sediment. 
Nissenbaum et al. (1972) identified the effects of methane production on the isotopic 
composition of EC0 2 . These studies were followed by numerous applications of 
isotopic measurements of SC0 2 to identify the processes that produce and utilize 
EC0 2 in sediments (e.g. LaZerte, 1981; Reeburgh, 1982; McCorkle et al., 1985; 

Herczeg, 1988; McCorkle and Emerson, 1988). 

Closed system box models were first used to try to reproduce the isotopic 
signatures that were measured in the field (Nissenbaum et al., 1972; Claypool and 
Kaplan, 1974; Reeburgh, 1982). Applying open system models to stable isotope 
studies of both freshwater and marine studies allowed for the inclusion of diffusion, 
bioturbation and sedimentation effects and yielded more information about the 
processes affecting these profiles such as carbonate dissolution (McCorkle et al., 1985; 
Alperin, 1988; McNichol et al., 1991), carbon input rates (McCorkle et al., 1985), 
and methane oxidation (Whelen et al., 1976; Reeburgh, 1982; Herczeg, 1988, 

Alperin, 1988). To date, diagenetic models of isotopic profiles have required 
assumptions about the isotopic composition of the UC0 2 produced during sulfate 
reduction and methanogenesis, as well as the isotopic composition of the remineralized 
organic carbon. 
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Typical SC0 2 6 13 C depth profiles from organic-rich marine sediments exhibit 
two major trends. The EC0 2 initially becomes depleted in 13 C relative to overlying 
seawater SC0 2 values (typically near zero per mil) and then, at depth, reverses this 
trend and becomes 13 C-enriched. The initial 13 C/ 12 C gradient has been attributed to 
sulfate reduction which produces SC0 2 with an isotopic composition similar to the 
organic matter that has been remineralized (typically -19 to -24 per mil in marine 
sediments (Fry and Sherr, 1984). Methanogenesis is thought to cause the gradient 
reversal by utilizing 13 C-depleted SCO ? and leaving behind 13 C-enriched SC0 2 
(Nissenbaum et al., 1972). To observe such profiles in typical coastal or deep sea 
sediments may require analysis of several to 100's of meters of sediment because of 
the slow rates of remineralization. The sediments of Cape Lookout Bight, North 
Carolina (Station A-l; Martens, 1976) were chosen for this study because these trends 
are fully established in the upper 40 cm of sediment due to the extremely high 
remineralization rates (Chapter 2; Boehme et al., in revision). 

The present study attempts to address some of the assumptions applied to 
isotopic studies of EC0 2 and to identify and quantify the controls on the isotopic 
composition of SC0 2 in a methanogenic sediment. Specifically, this study was 
undertaken to determine quantitatively the processes controlling the isotopic 
composition of SC0 2 in a methane producing marine sediment overlain by a typical 
sulfate reduction zone. By coupling field and laboratory studies, the parameters that 
had to be assumed in earlier models could be directly measured in order to estimate 
the relative importance of the dominant processes controlling the 6 13 C composition of 


53 



the £C0 2 . Field measurements of SC0 2 6 l3 C sedimentary profiles indicated a 
relationship between rates of the dominant processes and the resultant SC0 2 6 C 
profiles. Sediment incubations were used to determine rates of SC0 2 production and 
sulfate reduction as well as the 6 13 C of SC0 2 produced. From these measurements it 
was possible to test the hypothesis that the SC0 2 6 13 C porewater profiles from anoxic 
sediments resulted from a balance of SC0 2 produced by sulfate reduction and that 
produced by methanogenesis. 

Field Site 

Cape Lookout Bight, North Carolina is a partially enclosed marine basin 
located 110 km southwest of Cape Hatteras. Samples were collected at Station A-l, 
which has been actively studied for almost 20 years (Martens, 1976; Chanton et al., 
1983; Martens and Klump, 1984; and references therein). The dominant 
remineralization processes within the sediment are sulfate reduction and 
methanogenesis (Crill and Martens, 1983; 1986; Martens and Klump, 1984). 
Accumulation rates of 8-12 cm per year for the upper 40 cm of the sediment have 
been measured at Station A-l (Chanton et al, 1983; Canuel et al., 1990). Seasonal 
temperature variations of 20°C drive changes in remineralization rates and fluxes of 
the diagenetic products. Despite this seasonality, near-steady state conditions occur on 
yearly time scales. Good agreement between measured and modelled nutrient profiles 
suggest near steady state conditions may occur on shorter time scales as well (Klump 
and Martens, 1981; Martens and Klump, 1984; Haddad and Martens, 1990; Chapter 
2, Boehme et al., in revision; this study). 


METHODS 


The collection and analysis of field samples for monthly profiles of EC0 2 from 
Station A-l have been described previously (Chapter 2; Boehme et al., in revision). 

A 9.5 cm diameter sediment core was collected at A-l in June 1991 and sectioned in 
the lab at 2 cm intervals down to 20 cm to conduct an incubation experiment. 

Sediment temperature was 25 °C, in the field at the time of core collection. Sediments 
from each depth interval were homogenized and transferred into five 15 ml centrifuge 
tubes via syringes. The centrifuge tubes from each of the 10 depth intervals were 
capped and stored in a glove bag and kept under nitrogen at 25 °C for the duration of 
the incubation experiment. 

One sample from each of the ten depth intervals was processed immediately 
(T=0) and then every 40 to 75 hours by centrifugation (3500 rpm for 30 minutes). 

The extracted porewater was taken up in a syringe, and filtered through a 0.45 pm 
nylon filter (Micron Separations Inc.). One ml of porewater was injected into a 2.5 
ml serum bottle, capped and crimped for analysis of the concentration and 6 C of the 
EC0 2 . The EC0 2 samples were maintained frozen until the end of the experiment so 
that the five samples from each depth could be analyzed on the same day. One to 3 
ml of porewater was acidified and stored for Ca + + analyses. The remaining 
porewater was treated with a drop of concentrated ZnCl 2 to precipitate sulfide, 
filtered, and stored refrigerated for sulfate analysis. 

Porewater EC0 2 samples were analyzed on a GC-TCD adapted with a vacuum 
line to collect the C0 2 gas for isotopic analysis (Blair and Carter, 1992; Schaff et al.. 
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1993). The SC0 2 serum bottles were injected with 0.2 ml of 1M H 3 P0 4 and the 
resulting C0 2 was swept into the GC column containing Unibeads- Is or 2s silica gel 
(Alltech) with He as the carrier gas. Water was removed using a Nafion drying tube 
(Permapure Products). Sample concentrations were determined with a thermal 
conductivity detector. The C0 2 was trapped in an 1/8 inch o.d. stainless steel trap 
under liquid nitrogen, and transferred to 6 mm borosilicate glass tubing for storage. 
Isotopic analyses were done on a modified Finnigan Mat Delta E isotope ratio mass 
spectrometer (Hayes, 1983). All isotope values are given in the 6 C notation (See 
Eqn. 1). The precision and accuracy of the SC0 2 concentration measurements 

13 

determined from standards is 0.5 mM and 1.0 mM respectively, and for SC0 2 fi C 
measurements, 0.3 per mil and 0.6 per mil respectively. 

Porewater S0 4 = samples were processed by gravimetric analysis of the 
precipitated barium salt (Chanton, 1985; Chanton et al, 1987). The precision and 
accuracy of this technique was 0.6 mM and 1.0 mM, respectively. 

Porewater Ca + + samples were treated as described in Chapter 2. The Ca 
measure ments were analyzed on a Perkin-Elmer Atomic Absorption spectrometer in 
the Soil Science Laboratory (NCSU). 

RESULTS 

The concentration of SC0 2 increased and that of sulfate decreased as a function 
of time in the upper 10 cm during the incubation experiment, with the most rapid 
changes occurring in the upper 0-2 cm interval (Fig. 3.1). SC0 2 6 13 C became 
progressively depleted in * 3 C with time for the 0-8 cm intervals, whereas the 8-10 cm 
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interval became enriched in 13 C with time (Fig. 3.1). The rates of SC0 2 production 
and sulfate reduction were determined by linear fits to the concentration data for each 
of the upper 0-10 cm intervals (r 2 values are given in Table 3.1). The S0 4 - and 
EC0 2 concentrations, and the SC0 2 5 13 C values showed little change with time in the 
samples below 10 cm due to slow rates of remineralization and the relatively short 10- 
day incubation and therefore are not shown. The Ca + + concentrations show little 
change below the 0-2 cm interval (Fig. 3.5) for the upper 10 cm. The Ca" t ’ + T=0 
profile did not agree well with the porewater collected by squeezer on the same day 
(Chapter 2; Fig. 2.5). The poor agreement between the squeezer and centrifuged 
profiles may indicate that one of these techniques causes changes in the porewater 
concentrations, and therefore this data set was only used to show that there were not 
large changes in porewater Ca + + over the 10 day experiment. 

Using the EC0 2 concentration and 5 13 C values, 12 C and 13 C production rates 
were determined. Rate measurements of SC0 2 production, and S0 4 reduction for 
the upper 10 cm (Fig. 3.2a) as well as S 13 C0 2 and S 12 C0 2 production rates 
calculated using the isotopic signatures of the SC0 2 produced (Fig. 3.2b) were fit to 
an exponential curve of the form: 


R Z =R 0 exp (ai) 


( 10 ) 


where 

R z = rate of production or reduction of process (mM/hr), 
Rq = rate at the sediment-water interface (mM/hr), 
a = attenuation coefficient (cm’ 1 ), and 
z = depth into sediment (cm). 
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Values of Rq and a were solved iteratively from measured values of Rz with the 
Sigmaplot curve fitting routine (Jandel Scientific). Curve fitting parameters are given 
in Table 3.2. The EC0 2 concentration and 8 13 C values were used to estimate the 
6 l3 C signal of the EC0 2 produced for each depth interval using the following 
equation: 






'add 


(ID 


where 

C t = SC0 2 concentration at given time 
8 t = 6 13 C of EC0 2 at given time 

C = EC0 2 concentration at initial time zero for given depth 

5° = 8 13 C of EC0 2 at initial time zero for given depth 

c add = sco 2 concentration added between C 0 and C t 
5 add = 6 13 C of EC0 2 added between C Q and C t . 

The 6 aHa was solved for by curvefitting the EC0 2 concentration and EC0 2 8 C for 

each depth interval and the horizontal error bars represent the standard error for each 

of the values (Fig. 3.2). 

The advection-diffusion model (Berner, 1980) was used to convert the 
measured rates of sulfate reduction, EC0 2 production, E 12 C0 2 production and E 13 C0 2 
production to test whether these rates are representative of in siiU conditions. An 


analytical solution of the advection-diffusion equation 

X, D £C_ a dC_ R (12) 

dt s dz 2 

was found assuming dC/dt — 0 and the following boundary conditions: as C z — > °°, 
dC/dz = 0; for EC0 2 , C 0 = 2.11 mM, fi 13 C = 0.0 per mil; for sulfate, C 0 = 29.6. 
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Eqn. 12 becomes 


C Z =C 0 +— — — (1 -exp"), (13) 

arD s +a co 

where 

D = Diffusion coefficient (D s = <|) 2 D 0 ; <J)= 0.9 for upper 20 cm; 

Chanton et al., 1983; D s (HC0 3 ') = 0.0344 cm/hr; 

D s (S 0 4 = ) = 0.0312 cm/hr), 

C z = Concentration of species at given depth z 

a = attenuation coefficient (cm' 1 ) 

<u = sedimentation rate (10.6 cm/yr; Chanton et al., 1983). 

Diffusion coefficients for HC0 3 ' and S0 4 = were taken from Li and Gregory (1974) 
and corrected for sediment tortuosity as noted above (Ullman and Aller, 1982). The 
effects of changing porosity with depth have been shown to alter the results by less 
than 4% for Cape Lookout sediments and so porosity was assumed to be constant with 
depth (Klump and Martens, 1989). The comparison of measured and calculated 
concentrations and the EC0 2 6 13 C profiles are shown in Fig. 3.3. The steady state 
assumption appears to be justified given the agreement between measured and 
estimated EC0 2 profiles determined from the rate measurements (Klump and Martens, 
1981 ; 1989; this study). The estimated EC0 2 6 13 C profile shown in Fig. 3.3b, agrees 
with the measured values for the upper 0-12 cm. The two profiles diverge below this 
depth, because the lower boundary condition forces dCldz to 0 at depth. The 
modelled and measured S0 4 = concentration also agree well, verifying the sulfate 
reduction rate measurements. 

The depth-integrated rates of EC0 2 production and sulfate reduction as well as 
E l2 C0 2 and S 13 C0 2 were determined by integrating Eqn. 11 and multiplying by a 
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porosity term (<J>=0.9 cm^cm 3 ^ ) (Eqn. 12), to compare to previous depth 
integrated rate measurements: 

2 i 

The calc ulate d depth integrated rate of sulfate reduction of 21.3 moles-m -yr is m 
excellent agreement with a previously measured incubation experiment sulfate 
reduction rate of 21.0 moles-m^-yr' 1 (Klump and Martens 1989) and comparable to 
previously measured rates of sulfate reduction for June 1978 of 19.3 ± 3.4 moles-m 
2 -yr _1 and June 1979 of 24.3 ± 3.0 moles- m' 2 -yr‘ 1 (Martens and Klump, 1984). The 
depth integrated SC0 2 production of 34.6 moles-m^-yr' 1 also agrees well with a 
previous incubation experiment rate of 34.0 moles-m^-yr* 1 (Klump and Martens 
1989) and is comparable to EC0 2 fluxes of 38.9 + 0.8 moles-m 2 -yr 1 measured in 
June 1978 (Martens and Klump, 1984). 

DISCUSSION 

The EC0 2 6 13 C profiles measured at Cape Lookout are typical of anoxic 
marine sediments, with a zone of 13 C-depletion overlying a zone of * 3 C enriched 
EC0 2 . These profiles are thought to result from sulfate reduction producing EC0 2 
that is isotopically similar to the organic matter being oxidized in the upper zone and 
methanogenesis generating 13 C-enriched SC0 2 at depth. Additionally, 13 C-depletions 
have been attributed to methane oxidation adding * 3 C-depleted EC0 2 in the sulfate 
reduction zone, particularly near the transition between the sulfate reduction and 
methane production zones (Reeburgh, 1982; Whiticar and Faber, 1986). In studies 
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where methane oxidation has been shown to be important, the resultant EC0 2 may be 
distinguishable by the * 3 C-depleted nature and the concave-up sedimentary methane 
concentration profiles. At Cape Lookout, methane oxidation is not an obvious control 
on the SC0 2 6 13 C values, and appears to affect only the surficial sediments based on 
flux experiments and sedimentary CH 4 6 13 C profiles (See Chapter 2, Fig. 2.4; 

Boehme et al., in revision). The agreement between methane production rates and 
observed fluxes (Crill and Martens, 1986) is further evidence against significant 
methane oxidation at this site. 

The £C0 2 6 13 C profiles from Cape Lookout also exhibit seasonal changes 
(Chapter 2; Boehme et al., in revision) If sulfate reduction and methanogenesis are 
controlling the isotopic signature of the EC0 2 , then the seasonal trends seen in the 
isotope profiles should correlate to the changing rates of these processes. This can be 
tested directly using the measured monthly profiles of EC0 2 6 C. 

Further, if these processes are occurring in distinct zones, namely sulfate 
reducing bacteria out-competing the methanogenic consortia down to the depth at 
which sulfate is depleted, then the isotopic signal of EC0 2 being produced should 
reflect this zonation and should be expressed in the incubation experimental results. 

The measured isotopic signature for the EC0 2 produced in the sediment generally 
reflects this zonation (Fig. 3.2b). 

The use of exponential fits to the E 12 C0 2 and the E 13 C0 2 rate measurements 
results in an estimated 6^ 3 C produced that does not reflect the observed zonation in 
E 12 COo and E 13 C0 2 production (Fig 3.2b). This is not surprising considering that the 
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S 12 C0 2 and E* 3 C0 2 rate data may not be described by a simple exponential function 
b eca u se these rates result from the addition of two separate processes, sulfate 
reduction and methanogenesis, with different rates. The exponential fits, despite these 
problems reproduced the concentration and estimated isotope signatures of the EC0 2 
measured in the sediment (Fig. 3.3; see discussion below). A comparison of 
exponential fits to cubic spline fits resulted in little differences in the estimated EC0 2 
S l3 C profiles. The curve fits were used to determine depth integrated rates for 
modelling of individual processes. The depth integrated rates determined are 
consistent with previously measured rates for this site. 

Field Studies 

Seasonal changes in the HC0 2 6 13 C profiles from Cape Lookout (Chapter 2, 
Fig. 2.2; Boehme et al., in revision) appear to result from the changing rates of 
sulfate reduction and methanogenesis, the changing depths of sulfate depletion and 
possibly by changes in methanogenic pathways. This relationship can be seen 
qualitatively in the steepening concentration gradients and changing porewater SC0 2 
isotope profiles in Fig. 2.2. In an attempt to quantify this relationship, sulfate 
concentration gradients and the isotopic gradient for the upper 3 cm of the sediment 
are compared in Fig. 3.4. Increasing rates generally correlate with increasingly C- 
depleted gradients, however, this comparison does not include the effects of 
methanogenesis on the isotope profiles. Methane production is altering the isotopic 
signature of the SC0 2 flux. Further, for most of the year methane is produced 
predominantly via C0 2 reduction, however in July and August, some fraction of the 
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methane appears to be produced by acetate dissimilation (Crill and Martens, 1986; 

Blair et al., 1993). The smaller fractionation factor associated with acetate 
dissimilation may further alter the relationship between the sulfate concentration 
gradient and isotopic gradient for the months when acetate dissimilation is important. 

The comparison shown in Fig. 3.4 suggests a simple correlation, that SC0 2 
6 13 C values at the sediment surface are dependent on the concentration gradients of 
sulfate reduction. This may be a useful observation for understanding the range of 
isotopic signatures that are measured in other field sites, especially environments 
where methane production is not important. The data imply that if the isotopic 
signature of the two dominant processes controlling EC0 2 concentration in Cape 
Lookout sediments can be determined, then the overall SC0 2 isotopic signature 
should be predictable based on the changing rates of these two processes. 

Incubation Experiment 

The incubation experiment was performed to measure directly the isotopic 
signature of the EC0 2 being produced as a function of depth. Variation with depth in 
the isotopic composition of the EC0 2 produced for the upper 10 cm are shown in Fig. 
3.2b. For the upper 8 cm, the isotopic signature is relatively constant. The 8-10 cm 
interval is markedly 13 C-enriched with an isotopic signature of +2 ± 1.9 per mil. 

The two distinct zones are consistent with biogeochemical arguments that have 
theorized little overlap between the sulfate reduction and methanogenic zones. The 
isotopic signature of the EC0 2 in the upper 8 cm is enriched in 13 C relative to the 
metabolizable organic carbon being at Cape Lookout (-19 per mil; Chapter 2; Boehme 
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et al., in revision). This may indicate some methane production in the sulfate 
reduction zone. This hypothesis is consistent with previous measurements of methane 
production in the sulfate reduction zone (Crill and Martens, 1986). There are other 
possible causes for this enrichment of the SCO 2 that can be considered. 

One source of enriched SCO 2 to porewaters is the dissolution of calcium 
carbonate. Porewater concentrations of Ca + + were measured during the incubation 
experiment and suggest some dissolution in the upper 0-4 cm (Fig. 3.5). The previous 
Ca + + measurements for this site also show an increase in the upper two cm indicating 
dissolution (Chapter 2, Fig. 2.5; Boehme et al., in revision). The increases, 
however, are not large enough to account for an offset of -5 per mil (-19 per mil 
organic carbon to the -14 per mil average isotope value for the upper 8 cm) assuming 
an isotopic composition of the CaCOj of -0. 1 per mil (the average isotopic 
composition of the particulate inorganic carbon at Cape Lookout). Further, the Ca 
profiles measured during the incubation experiment do not indicate dissolution below 
the 0-2 cm interval over the 10 day incubation experiment (Fig. 3.5). Thus, the 
isotopic composition of the ECO 2 produced should therefore not be significantly 
affected by calcium carbonate dissolution. This conclusion could be erroneous if 
dissolution was coupled with precipitation in an exchange reaction. Exchange 
reactions have been hypothesized to account for differences in modelled and measured 
SC0 2 S l3 C profiles from turbidites (McArthur, 1989) but have typically been ignored 
be ca use they are considered to be too slow to significantly alter the isotopic 
composition of rapidly depositing sediments (Nissenbaum et al., 1972). The isotopic 
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composition of the particulate inorganic carbon is relatively constant with depth 
(Boehme et al., in revision; Boehme, 1989), however, the possibility of isotopic 
exchange reactions during dissolution and precipitation of CaC 03 has not been 
adequately addressed. 

Another possible cause of the relatively ^C-enriched SC0 2 in the sulfate 
reduction zone is the preferential remineralization of a ^C-ennched organic carbon 
fraction during sulfate reduction. A long-term incubation study using sediments from 
Cape Lookout estimated the 6^C of the organic matter remineralized during sulfate 
reduction to be -15.6 per mil (Alperin et al., 1992). If the sulfate reducing bacteria 
do not significantly fractionate the organic matter during remineralization, then the 
SC0 2 should be similar to the remineralized organic carbon. Preferential 
remineralization is not consistent with the particulate organic carbon (POC) 6 C 
profiles for this site that indicate that the 6 value of the organic carbon does not 
change with depth (Haddad, 1989; Blair and Carter, 1992; Chapter 2; Boehme et al., 
in revision). If sulfate reduction, the dominant remineralization process in these 
sediments, is preferentially removing organic carbon with a different isotopic signature 
than the bulk organic carbon, then this should cause be evident in the POC 5 C 
profiles. 

1 3 

If methane production is occurring within the sulfate reduction zone, C- 
enriched EC0 2 would be added to the pool of SC0 2 produced by sulfate reduction. 
This is consistent with incubation measurements of methanogenic rates measured by 
Crill and Martens (1983; 1986) in which low rates of methanogenesis via C0 2 
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reduction were measured within the sulfate reduction zone using incubation 
experiments and tracers, especially in the summer months when sulfate is rapidly 
depleted by 8 to 10 cm. As will be shown later, these rates of sulfate reduction and 
methanogenesis can be used in conjunction with an estimate of the isotopic signature 
of the ECO 2 produced from these processes to estimate a ECO 2 value produced 
in the s ulfate reduction zone. Studies have shown that methane production via non- 
competitive substrates such as methanol and methylamines can also result in methane 
production (and presumably EC0 2 production as well) in salt marsh sediments (King, 
1984 ; King et al., 1985 ; Oremland et al., 1982 ; 1993 ). 

The successful use of the rate measurements to estimate EC0 2 concentration 
and 6 13 C profiles indicates that a simple two component mixing model can be used to 
estimate the isotopic signature of the EC0 2 produced from an anoxic sediment like 
Cape Lookout. The mixing model assumes that the isotopic signature of the ECO 2 
results from the mixing of two sources— oxidative source (sulfate reduction) and 
methanogenic. This hypothesis can be described mathematically as a mass balance of 
the processes producing EC0 2 in the sediment, 

(15) 

EC0 2 total production rate (mM/hr), 
isotopic signature of EC0 2 produced, 

EC0 2 production from sulfate reduction rate (mM/hr), 
isotopic signature of EC0 2 from SR 
EC0 2 production from methanogenesis rate (mM/hr), 
isotopic signature of ECO 2 from methanogenesis. 


where: 

Rtot ~ 

^ tot = 

R(EC0 2 ) sr = 

5 SR = 

R(EC0 2 ) m = 

6m 
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The sources of the 13 C-enriched EC0 2 produced in the sulfate reduction zone cannot 
be resolved by this study, but the measurement of this isotope signature allows us to 
model the processes producing these signals in an effort to test some of the hypotheses 
given here. 

The SC0 2 produced and its isotopic signature were used to determine the 12 C 
and 13 C production rate profiles individually. These profiles were curve fit and 
extrapolated to estimate the isotopic signature of the SC0 2 added to the surface 
sediments due to sulfate reduction, -19.2 per mil. This value is assumed to be the 
best guess value for the isotopic composition of the SC0 2 produced from sulfate 
reduction alone (fi SR ) because the surface sediments should be the least affected by 
methanogenesis. 

The sulfate reduction rate, R SR , can be used to determine the depth integrated 
rate of SC0 2 produced from sulfate reduction, 

R(HCO^ sr =R sr z, 

where z — 1.78, the ratio of SC0 2 produced to sulfate reduced at the sediment-water 
interface. The stoichiometric coefficient, z, is assumed to be constant downcore. 

This assumption is supported by degradation studies of algae where the calculated 
composition of the refractory material (calculated as percent proteins, lipids and 
carbohydrates) was similar to the original composition of the algae (Foree and 
McCarty 1970) suggesting that no one algal component was preferentially utilized 
producing a different ratio of SC0 2 produced to algae degraded. Previous estimates 
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of t for Cape Lookout surface sediments, of 1.7 ±0.1 (Alperin et al, 1992) and 1.9 
(average of upper 6 cm; Klump and Martens, 1989) supports the x value used in this 
study. Factors controlling t are discussed in further detail in a later section. 

The net production of EC0 2 during methanogenesis, R M is calculated by 
subtracting R(SC0 2 ) sr from R tot . The depth integrated 5 value of SC0 2 produced, 
6 tot was estimated from the isotopic signature of the depth integrated 2 C0 2 and 
S 13 C0 2 production rates, -13.1 per mil. Solving Eqn. 15 for 6 M gives +44.2 per mil 
for the isotopic composition of the SC0 2 produced during methanogenesis. Using the 
values for SC0 2 production from sulfate reduction (-19.2 per mil) and methanogenesis 
(+44.2 per mil) and the estimated production rates, a SC0 2 fi l3 C mixing curve was 
calculated and compared to the measured isotopic signature of the SC0 2 produced 
(Fig. 3.2b). The mixing curve misses the zonation of sedimentary processes, similar 
to the exponential curve fit results, but the mixing curve is derived from the 
exponential fits and so the similarity is not unexpected. As noted previously, the 
estimated concentration from the rate curve fits agrees well with the data. The 
differences in the curve fits and the measured data are probably not strongly reflected 
in the estimated isotope profiles because the production rates are low at depth. 

Based on mass balance arguments, it should be possible to estimate the isotopic 
composition of the CH 4 produced in the incubation experiment using the following 
equation, 

+ ^ remix Cff 4 ^ 
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where 


R CH4 = rate of methanogenesis (mM/hr), 

5 CH4 = isotopic signature of methane produced, 

5 remi n = (-18.9 per mil; Boehme et al., in revision). 

To estimate the rate of methane production, the rate of ECO 2 produced during 

methanogenesis can be used in the following stoichiometric relationship. 

(R R{1,CO^ m (18) 

{ CB * ] (Ti T‘ 

The x- 1 parameter was determined based on oxidation state arguments and will be 
discussed later. Solving Eqn. 17, the isotopic composition of the methane produced in 
the sediments is -65.9 per mil. This value is very similar to the measured isotopic 
signature of methane bubbles collected at the same site in June 1983 and 1984 (-64.3 
± 0.7; Martens et al., 1986). The mass balance equations accurately estimate the 
isotopic composition of the methane being produced in this sediment. This is further 
verification that the mass balance achieved using measured rates of SC0 2 production 
and sulfate reduction adequately represents controlling diagenetic processes. 

As noted earlier, the calculated 6^C for EC0 2 produced via sulfate reduction 
(-19 per mil) and methanogenesis (+44 per mil) can be used to estimate the isotopic 
signature of the SC0 2 produced in the sulfate reduction zone using previously 
measured rates from Crill and Martens (1983). The rates of sulfate reduction and 
methane production from the summer months were averaged for the upper 1-4 cm 
and 6-1 lcm (the depths above where sulfate was depleted) and converted to £C0 2 
production rates using x and x-1. The resultant rates were used with the calculated 
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$ 13 C for EC0 2 produced via sulfate reduction (-19 per mil) and methanogenesis (+44 
per mil) to estimate a 6 13 C for SC0 2 produced of -18 per mil for the upper 4 cm and 
-13 per mil for the 6- 11 cm interval. These calculations suggest that methane 
production in the sulfate reduction zone can significantly affect the 6 C of the EC0 2 
produced and may be responsible for the 13 C-enriched EC0 2 seen in the sulfate 
reduction zone. 

The determination of the isotopic composition of the methane produced in the 
sediment allows us to calculate the fractionation factor for the production of methane, 

(19) 

5 ,3 W I® 3 ’ 

where 6 l3 C( CH4 ) is -65.9 per mil as calculated above and 6C0 2 is the isotopic 
composition of the C0 2 measured in the methanogenic zone. Using a pH of 6.95 for 
A-l sediments (Chanton, 1985; N.E. Blair et al., 1993), the relative contributions of 
HCO3', C0 2 , and C0 3 = can be determined (Stumm and Morgan, 1981). At 8-10 
cm, the measured SC0 2 6 13 C is -7.0 per mil (Fig. 3.3). To determine the isotopic 
composition of the C0 2 species, the following equations (Deines et al., 1974; 
Friedman and O'Neil, 1977; Blair et al., 1993) were solved simultaneously: 

6 13 C(SC0 2 )=d6 13 C(C0 2 )+e6 13 C(HC0 3 ')t/3 13 C(C03’) (20) 

HCO 3 lOM 13 C(ffCQ 3 ~) (2i) 

“ C0 2 ’ 1( J 3 +6 13 C(C0 2 ) 
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( 22 ) 


5 13 C(ffCO 3 >S 13 C(C0 3 ’) 

The variables d, e, and f represent the fractions of the dissolved SC0 2 components. 

The 6 13 C of HC0 3 ' and C0 3 = are assumed to be the same (Eqn. 22). The 
equilibrium fractionation factor for the HC0 3 ’ and C0 2 equilibrioum is given by: 

ln a (^^)=(^^) -0.0241 (23) 

k co 2 T 

Temperature is in Kelvin (Deines et al., 1974) Solving Eqns. 20, 21, 22 and 23 gives 
a 6 13 C value for C0 2 of -14.0 per mil. Substituting this value into Eqn. 19 gives an 
a for methane production of 1.055. This fractionation factor is remarkably similar to 
an independent estimate for this site, based on a mechanistic model (1.056 at 25 C; 
Blair et al., 1993) and is within the range determined in culture studies (1.03 to 1.06; 
Games et al., 1976; Fuchs et al., 1979; Balabane et al., 1987; Belyaev et al.,1983) 
and estimated from measured profiles of SC0 2 and CH 4 (1.05 to 1.09; Whiticar et al., 
1986). The agreement is further evidence to support the model results. This is the 
first in situ fractionation factor for methane production from C0 2 reduction 
determined using directly measured rates and the associated isotopic signatures. 

Model Sensitivity to Reaction Stoichiometry 

The isotopic composition of the SC0 2 produced during sulfate reduction and 
methanogenesis, and the isotopic composition of the methane that were estimated using 
the mass balance approach are sensitive to the stoichiometries of the sulfate reduction 
and methane production reactions. The 6 C jj 4 is especially sensitive to the 
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stoichiometry as is seen in Fig. 3.6. These stoichiometries are dependent on the 
apparent oxidation state of the organic carbon that is remineralized. For the 
generalized formula C x HyO z , the apparent oxidation state of the carbon is defined as 


for the reaction 


OS=^ 1 for-4iOSi4. 

x 


(24) 


C x H y 0 z +wS0^uC0 2 +S ‘ +tf 2 0, (25) 

where 

x = (26) 

w 4-OS 

Similarly for the methanogenic reaction, 

2C x H y 0-uC0 2 +vCH A (27) 

the ratio of C0 2 produced to CH 4 produced (u/v) is given by 

u 4+OS 

v~ 4-OS~ 

The relationship between t-1, the ratio of CH 4 to EC0 2 produced, and oxidation state 
has been demonstrated in culture experiments (Tarvin and Buswell, 1934; Fig. 3.7). 

Solving Eqn. 28 for oxidation state using t = 1.78 from the model, gives an 
apparent oxidation state of the organic matter of -0.49. Using the t determined by 



Alperin et al. (1992) of 1.7 ± 0.1 gives an OS values of -0.72 ± 0.28 and using the 
x determined by Klump and Martens (1989) of 1.9 gives an OS value of -0.21. 

An independent estimate of oxidation state of the remineralized organic matter 
can be determined from the composition of the remineraiized organic matter. At Cape 
Lookout, 64 ± 17% of the metabolized organic carbon has been identified as 
carbohydrates, lipids and amino acids, in a ratio of l-Otl.0: 1.9. Using the equation, 

os= (3a+2b+2z~y) (29) 

x 

to describe the oxidation state of amino acids with the formula C x H y N a O z Sjj and the 
amino acid distributions identified at Cape Lookout sediments (Burdige and Martens, 
1988), the average oxidation state of the amino acid carbon is estimated to be 0.07. 
The carbohydrate carbon is assumed to have an average oxidation state of zero. The 
C l6 fatty acid was used as a representative lipid for this system (Haddad, 1989) with 
an oxidation state of -1.75. Given the relative contribution of these fractions to the 
identified pool, the apparent oxidation state of the identified metabolizable pool is 
-0.41, in good agreement with the modeled value of -0.49. The OS determined from 
the organic matter can be applied to Eqns. 26 and 28 to estimate a x of 1.81, in 
excellent agreement with the 1.78 determined from this study. McCorkle and 
Emerson (1988) calculated C0 2 /0 2 ratios for a variety of oxic and suboxic sediments. 
Assuming that oxidation state is the only factor controlling this ratio, their value of - 
0.54 is consistent with the t determined for Cape Lookout. Gujer and Zehnder 
(1983) also noted a relationship between oxidation state, substrate, and C0 2 to CH^ 
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ratios for anaerobic digesters. Based on their correlations, the substrates that correlate 
to the oxidation state estimated for Cape Lookout are algae, bacteria, carbohydrates 
and proteins. The agreement of these various approaches to determining the OS and 
stoichiometries of sulfate reduction and methanogenesis is further evidence that despite 
the sensitivity of the model to these parameters (Fig. 3.6), our model is accurately 
describing the processes and their isotopic signatures in this system. 

Implications for Control of Marine CH 4 6 C 

Marine CH 4 6 13 C values exhibit a wide range of values (-110 to -60 per mil; 
Whiticar et al., 1986). Possible sources of this variation must be due to differences in 
the relative rates of methanogenic pathways (predominantly C0 2 -reduction and acetate 
dissimilation), differences in fractionation factors, and the isotopic composition of the 
precursors. The mass balance model presented here indicates that the 6CH 4 is also 
dependent on the relative rates of non-methanogenic oxidative processes (sulfate 
reduction in this study) and methanogenesis. Further, the model's sensitivity to 
oxidation state suggests that oxidation state may be important as well. 

CONCLUSIONS 

The EC0 2 fi 13 C profiles from Cape Lookout have been shown to be dependent 

on the rates of the remineralization processes, similar to the results of McCorkle and 

13 

Emerson (1988). In the methanogenic sediments of Cape Lookout, the EC0 2 6 C 
primarily reflects changes in sulfate reduction rate and methanogenesis. The 
sedimentary profiles are best described by the mixing of EC0 2 from two processes: 
sulfate reduction, producing EC0 2 with an isotopic signature of -19.2 per mil, and 



methanogenesis, resulting in EC0 2 with an isotopic signature of +44.2 per mil. The 
mass balance calculations generated a reasonable isotope value for methane and a 
reasonable fractionation factor for methane production. However the calculation is 

very sensitive to the reaction stoichiometry. 

The mass balance approach used suggests that the isotopic signature of both the 
SC0 2 and methane are sensitive to the ratios of EC0 2 /SR and SC0 2 /CH 4 that are 
ultimately dependent on the oxidation state of the organic matter being remineralized. 
The oxidation state determined from the incubation experiment agreed with a 
calculated OS based on the identified fraction of the organic carbon remineralized. 

We hypothesize that the ratio of sulfate reduction to methanogenesis may be an 
important contribution to the range of CH 4 isotope values observed in the marine 
environment. This may ultimately be useful in interpreting SC0 2 6 13 C profiles in 
other environments as well as the isotopic signature of diagenetic carbonates. 
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Fig. 3.1 SC0 2 , sulfate and SC0 2 6 13 C data for the five depth intervals. The 

EC0 2 , E 12 C0 2 , E 13 C0 2 and sulfate concentrations were fit to a line to 
determine rates for each depth interval. 
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Fig. 3.2 (A) Measured rates (symbols) and curvefits to rates (lines) of SC0 2 

production and sulfate reduction. (B) Isotopic signature of the SC0 2 
produced in the upper 10 cm (symbols) determined using Eqn. 11. 

The solid line is the 6 13 C of the ratio of the curvefits of the 
individual 2 12 C0 2 and 2 13 C0 2 production profiles. The dotted line 

( ) is a mixing curve based on the calcualted rates of sulfate 

reduction and methanogenesis and the associated end member isotopic 

signatures. The dashed line ( ) is the ratio of depth-dependent 

production rates of 2 *^C0 2 and 2 13 C0 2 solved for using the measured 
SC0 2 concentration and fi* 3 C profile and the diagenetic equation 
(Eqn. 12). 
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Fig. 3.3 The measured concentrations of EC0 2 and sulfate (symbols) and the 
model estimates determined using Eqn. 12 and the rate data shown in 
Fig. 3.2 (solid line) (A). The measured isotopic signature of the EC0 2 
(symbols) and the modelled estimates determined using Eqn. 12 and the 
individual rates of S^C0 2 and S^C0 2 production (solid line) (B). 
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Fig. 3.4 The sulfate gradient versus the EC0 2 6 13 C gradient for the upper 3 cm 
of sediment from porewater profiles collected from 1986 to 1991 at 
Cape Lookout. The sulfate and SC0 2 6 13 C gradients were based on 
linear fits of concentration over the upper 3 cm (one cm intervals). 
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Fig. 3.5 Porewater Calcium concentration versus time for the upper 10 cm from 
the incubation experiment. 


O 1, 




(w m ) ++ «0 


85 


Time (days) 








Fig. 3.6 The relationship between the oxidation state of the organic carbon, the 
ratio of SC0 2 produced to sulfate reduced (t) and the modelled isotopic 
signature of the methane produced. 
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Fig. 3.7 The relationship of the relative rates of C0 2 and CH 4 production ( t-1) 
and the apparent oxidation state (OS) of the fermented material. Data is 
from Tarvin and Buswell (1934). Curve is described by Eqn. 18. 
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Table 3.1 


13 

Regression coefficients (r^) for the S0 4 , SCO 2 , and SCO 2 6 C 
incubation data 


Depth 

so 4 = 

O 

u 

« 

am 

0-2 cm 

0.92 

0.99 

0.99 

2-4 cm 

0.94 

0.98 

0.98 

4-6 cm 

0.75 


0.96 

6-8 cm 

0.85 

0.89 

0.89 

8-10 cm 

0.56 

0.84 

0.84 


QO 
















Table 3.2 


Curve fit parameters, estimated depth integrated rates (10 cm), and 
associated isotopic compositions from tube incubation experiment. 


Process 

Ro 

a 

"Depth 

Integrated 

Rate 

fi l3 C 

Total EC0 2 
Produced 

0.103 

0.2037 

0.439 

-13.1 

EC0 2 Produced 
from Sulfate 
Reduction 

0.103 

0.2341 

0.397 

-19.2 

EC0 2 Produced 
from 

Methanogenesis 



0.0420 

+44.2 

Methane Produced 

— 

— 

0.0539 

-65.9 

Sulfate Reduction 

0.058 

0.2341 

0.223 

i 

1 


*Units of mmol*cm‘ 2 *hr' 1 
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